Configuration-Space Location of the Entanglement between Two Subsystems 



O 
O 

(N 

a 

o 

m 

(N 
> 
00 

l> 

O 

o 
o 

Oh! 

> 



H.-C. LiiE and A.J. FisheiQ 

UCL Department of Physics and Astronomy and London Centre for Nanotechnology, 
University College London, Gower Street, London WCIE 6BT, U.K. 

In this paper we address the question: where in configuration space is the entanglement be- 
tween two particles located? We present a thought-experiment, equally applicable to discrete or 
continuous-variable systems, in which one or both parties makes a preliminary measurement of the 
state with only enough resolution to determine whether or not the particle resides in a chosen region, 
before attempting to make use of the entanglement. We argue that this provides an operational 
answer to the question of how much entanglement was originally located within the chosen region. 
We illustrate the approach in a spin system, and also in a pair of coupled harmonic oscillators. Our 
approach is particularly simple to implement for pure states, since in this case the sub-ensemble in 
which the system is definitely located in the restricted region after the measurement is also pure, and 
hence its entanglement can be simply characterised by the entropy of the reduced density operators. 
For our spin example we present results showing how the entanglement varies as a function of the 
parameters of the initial state; for the continuous case, we find also how it depends on the location 
and size of the chosen regions. Hence we show that the distribution of entanglement is very different 
from the distribution of the classical correlations. 

PACS numbers: 03.67.Mn,03.65.Ud,42.50.Dv 



INTRODUCTION 



Studying the entanglement properties of a number of 
spatially extended many-body systems including spin 
chains, coupled fermions, and harmonic oscillators [J, 
i, i, S fi i, 0, B i, has both given informa- 
tion on the potential uses of these systems in quantum 
information processing, and yielded insight into their 
fundamental properties. Quantum entanglement is a 
measure of essentially quantum correlations, and many 
interacting systems possess an entangled ground state 

In this paper we address the question: where in con- 
figuration space is the entanglement between two parti- 
cles located? We pose the question using the language 
of spatial entanglement, which plays a significant role 
in many physical realizations of QIP (quantum informa- 
tion processing). However our results are easily recast 
in terms of other types of entanglement. Specifically, we 
investigate the location dependence of the ground-state 
entanglement between two interacting subsystems. We 
choose a pair of coupled harmonic oscillators as an ex- 
ample, since this is a system for which many exact results 
are available [1, [lB|- We assign one oscillator to each of 
the two communicating parties Alice and Bob, but per- 
form a thought experiment in which one or both of them 
first measure the system in configuration space, with just 
enough precision to localise it in some chosen region, and 
thereafter are restricted to operations only within that 
region. This restriction corresponds to a particular type 
of projective filtering in configuration space. We ask how 



this restriction affects the spatial entanglement available 
to them for other purposes — for example, for teleporting 
additional qubits between them. This should be distin- 
guished from the approach taken recently by Cavalcanti 
et. al. , who explored the effect of a finite-resolution 
spatial measurement on the spin entanglement of a sys- 
tem of noninteracting fermions and also a photonic in- 
terferometer. Our research also contrasts with previous 
studies (isl . [l9l . [20| of the entanglement of a finite region 
of space with the rest of the system. 

In a previous paper [2l| we investigated the limiting 
case where the size of the preliminary-measurement re- 
gion is very small, and showed that a smooth two-mode 
continuous-variable state can be approximated by a pair 
of qubits and its entanglement fully characterized, even 
for mixed states, by either concurrence density or nega- 
tivity density; here we shall focus on studying the varia- 
tions of the entanglement properties with the size of the 
region. For the present we assume that the two particles 
are distinguishable; the effects of indistinguishability on 
the phenomena discussed here are a subject for further 
work. We argue that the shared entanglement remaining 
to Alice and Bob provides a natural measure of where in 
configuration space the entanglement was originally lo- 
cated. We show that the distribution of entanglement is 
very different from that of the classical correlations. 



II. THEORY 



Restricting configuration space by Von 
Neumann measurements 



Let the configuration space of the whole system be de- 
* Electronic address: lho.lin@ucl.ac":iIkl scribed by the Coordinates qA and qs, where qa describes 

tElectronic address: ' andrew.fisher@ucl.ac.ukl Alice's particle and qs describes Bob's. We will initially 
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present the case in which only Alice makes a prepara- 
tory measurement on her system; suppose she has access 
to some restricted portion A of the configuration space 
of "her" particle, whose coordinate is qa- If she mea- 
sures her system with just enough accuracy to determine 
whether it is in region A or not, but no more, the effect is 
to localise the wavefunction either inside, or outside, the 
chosen region. The restriction to lying inside the region 
corresponds to the projector 



Ea^ \(1a) (qaI dqA (81 lother, 

Ja 



(1) 



where lothor is the identity operation for all the other 
particles (assumed distinguishable) in the system. 



1. The discarding ensemble 

Suppose A is of finite extent, and Alice measures the 
position of her particle with just enough accuracy to de- 
termine whether it is in A or not. If so, she keeps the 
state for further use; if not, she discards it (and tells Bob 
she has done so). Then the density matrix appropriate 
to the ensemble of retained systems is 



Pd,a 



EapEa 



(2) 



Tt{Eap) 

0AiqA)piqA,qothcr; q'a, gpthcr)^^ (g^) 
/ Ia P(lA, gothor; qA, gother)dqAd(7othcr ' 

where 9a is a generahzed Heavyside function defined so 
that 



OA{q) 



1 ifqeA 
otherwise 



(3) 



The subscript D refers to the discarding of the unwanted 
states; we refer to this density matrix as describing the 
"discarding ensemble". Note that, if the original p was a 
pure state Itfj) (V'l, then the post-selected density matrix 
is also pure: 



Pd,a 



Ea W (V'I Ea 
OjIEaW 



(4) 



In particular this means that even though the system has 
continuous variables and is therefore infinite-dimensional, 
its entanglement £d,a is easily calculated through the 
von Neumann entropy S'(p^'^-') of the reduced density ma- 
trix p(^) = Trsp. 



2. The non-discarding ensemble 

On the other hand if Alice chooses not to discard the 
system when she fails to detect a particle in region A, 
the appropriate density matrix is 



where the subscript ND refers to "non-discarding" and 
the complementary projector Ea' is defined as 



Ea' = 1- Ea 



qA^A 



I^a) (gyll dgyl (g) iothcr- (6) 



Eq. (lU describes a mixed state. It differs from the orig- 
inal density matrix p in that off-diagonal elements of p 
connecting qA &^A and qa ^ A have been set to zero. 

Let PA = T^t^[EapEa] be the probability of finding Al- 
ice's particle in A. Since the first and second components 
of pnd can be distinguished by Alice and Bob using local 
operations and classical communication (LOCC), they 
can teleport pa^d.a + (1 ^ Pa)£d.A' qubits on average 
between them. Hence the distallable entanglement (and 
therefore also the entanglement of formation) of pnd is 
not less than pa£d,a + (1 — Pa)£d,A'- On the other 
hand, Eq. ([5]) also constitutes a valid decomposition of 
the non-discarding density matrix pnd into orthogonal 
pure states; it follows that the entanglement of formation 
End is not greater than the average entanglement of this 
decomposition: End < Pa^d.a + (1 — Pa)Ed,A' ■ The 
only way these two observations can be consistent is if 



End = paEd,a + (1 - Pa)Ed,A'- 



(7) 



If all the operators available to Alice have support only 
in region A (i.e. if she can neither measure her particle's 
properties, or manipulate it in any way, except when it 
is in A) then the component projected by Ea' is "out of 
reach", and the second component Ea' pEa' of the state 
Pnd is functionally equivalent to a separable state as far 
as any operation that Alice and Bob can perform is con- 
cerned. It does not possess any entanglement properties 
that are useful to Alice and Bob. In that case, Eq. ([7]) 
reduces to End — PaEd.a- 



3. Precise measurements of 



If, on the other hand, Alice measures the position ac- 
curately, but again keeps only those occasions when the 
results lie within A, the discarding ensemble's density 
matrix is 



PP 



Ia EqpEqdq 
J^TT{E,p)dq' 



(8) 



Pnd = EapEa + Ea'pEa' 



(5) 



where the subscript P refers to measuring precisely and 
Eq is the projector corresponding to measuring Alice's 
particle A precisely at position q: 

Eq - S{qA - g). (9) 

Eq. ([H]) describes a density matrix that is diagonal in qA', 
it is a mixed state even if all the measurements where the 
particle is not found in A are discarded. Furthermore, 
unless there are some additional degrees of freedom of 
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particle qa which are not measured, the overall density 
matrix can be written as an incoherent sum of product 
states: 



where 



pp = EqpEqdq 
I A 



Iqa) {qA\dqA, 



A 



(10) 



where Iqa) is a state in which particle A is located ex- 
actly at QA and particle B is in some arbitrary state, pp 
therefore contains no remaining entanglement with Bob's 
particle B. 

Note that in the limit of very small measurement re- 
gions, the distinction between precise and imprecise mea- 
surements disappears. The case of vanishingly small re- 
gions was analysed in a previous paper [2lj , where it was 
shown that a well-defined concurrence density exists, and 
so the concurrence after the measurement is directly pro- 
portional to the region size. 



4- Measurements by both parties 

Exactly analogous formulae can be written down for 
the cases where Bob makes a preliminary measurement 
on his particle, or both partners make a measurement. 
In the case where both parties make a preliminary mea- 
surement, the reduced density matrix of Alice's system 
that is used to calculate the entanglement will naturally 
depend also on the measurement performed by Bob. 



PAB = Tr[£'Hi?^p] 



(15) 



is the probability of finding Alice's part of the system in 
A and Bob's part in B, and 



Pd,ab = 



PAB 



(16) 



is the density matrix in the discarding ensemble after 
this measurement result has been obtained. Combining 
Eq. lfT3|) and Eq. (fT4|) we obtain the following inequal- 
ity for the average of the entanglement in the discarding 
ensemble over all the partitions: 



^PAB£DipD,AB) <£{P)- 
AB 



B. Spin systems 



(17) 



We can make an exactly analogous theory for the case 
where Alice and Bob share a system defined on some 
other state space, for example a spin system — perhaps 
more familiar in quantum information theory. We simply 
replace the projection operator Ea by one defined in spin 
space; for example, Ea might project onto states with a 
specified spin component in a given direction. The rest 
of the theory is as outlined above. 



III. THE PHYSICAL SYSTEMS 



An inequality for the discarding entanglement 



Spins 



Suppose Alice and Bob divide their configuration 
spaces into a set of segments A and B respectively, and 
each make a measurement determining in which segment 
the system is located. In the nondiscarding ensemble, 
Eq. ^ generalizes to 



Pnd 



Y^EbEapEaEb, 



AB 



where 



(11) 



(12) 



However, this corresponds to a local operation performed 
by Alice and Bob. Their shared entanglement is non- 
increasing under this operation; therefore. 



^ip) > £{pnd)- 



(13) 



But, by a straightforward extension of the argument 
given above. 



£{pnd) = '^Pab£d{pd.,ab)i 



(14) 



Suppose that both Alice and Bob each possess two 
spins; the first spins belonging to each of them are en- 
tangled, as are the second spins, and the overall state 
of the system is a product of the state of the two pairs. 
For example, we could write 



1^/)) = (cos6'i ItAiTsi, 

(8)(C0S6'2 |Ta2Ts2) 



sin^i liAiisi)) 

sin02|U2iB2)); (18) 

the state is pure so entanglement between Alice's and 
Bob's subsystems is well quantified by the von Neumann 
entropy of the reduced density matrix. Suppose also that 
Alice and Bob can only handle systems if the total spins 
Ms available to each party are such that Ms = 0; per- 
haps the parts of the state with non-zero moment are lost 
because of the presence of large fiuctuating fields in the 
environment. In the discarding ensemble defined by this 
restriction, the state becomes 

1 



AB 



^/l — cos 29i cos 26*2 

(cos 6*1 ITaiTbi) sin 6*2 |iA2iB2) 

-|-sin6li lUiisi) 0036*2 |Ta2Tb2)), (19) 
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again this restricted state is pure but entanglement 
should become quite different. 

This type of measurement is familar in other 
contexts — for example entanglement distillation and con- 
centration [13, HI- 

Entanglement of a bi-partite mixed spin state can also 
be easily quantified by using negativity M{p) instead as 
the entanglement measure. We define negativity as the 
sum of the magnitudes of the negative eigenvalues of 
the partially transposed density matrix fF" , 



Af{p) 



iS.t. Ai<0 



Consider the mixed state defined by 

l-F ■ 



p = — — w) m 
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-1, 



(20) 



(21) 



where is as defined in Eq. (flSl) (in contrast to the 
definition of Werner states, this is not a maximally en- 
tangled state) and F S [1/16, 1]. Note that when F = 1^ 
the state becomes pure. Again an example of discarding 
ensembles can be obtained by projecting the state l(2T|) 
onto Ms = subspace and renormalizing accordingly. 



B. Harmonic oscillators 

The density matrix of a Gaussian state can be written 
in the coordinate representation [23| as 

(<?|pk') = p{t,<i) 

-\{q-qrm<i-i') 

+ '-{q-qm<l + ci% (22) 

where Ci is a normalization constant, and where L, M 
and K are real iV-dimensional matrices with L and M 
symmetric, while K is arbitrary. These matrices are re- 
lated to the covariance matrix 7 by 



2 ' 



loypL \/io 

-K 1 j i i(L + M)-i j I -K 1 



(23) 

We note that for a pure state, M = and K is symmetric. 

Consider a harmonic system with a Hamiltonian (tak- 
ing h= 1) 

H-R i In /{2m) (^4) 

where the vector R of quadrature operators is given by 
the positions Rj = Xj and conjugate momenta Rn+j = 
P], for I < j < N„ the N x N matrix V contains the 
coupling coefficients, and uj is the natural frequency of 



uncoupled oscillators. For a translationally invariant sys- 
tem the potential matrix elements depend only on the 
difference between the indices: Vj^fc — w for 

^ < j,k < N. The covariance matrix of the ground state 
is then Q| 
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-( ® muj-/p) 

2 mu! 



2 muj 



1/2 



® mtjVi/2). 



(25) 



Since the Hamiltonian given in Eq. I|24p has no coupling 
between position and momentum variables, 7 is block 
diagonal and hence K = 0. Furthermore if there are 
only nearest-neighbour interactions, with a Hooke's-law 
spring constant K, the interaction strength is character- 
ized by the single dimensionless parameter 

2K , , 

a = ^. 26 

For the two-oscillator ground state we therefore have 
only one non-zero matrix: 



1 + yi + 4a 1 - Vl + 4a 
1 - Vl + 4q; 1 + Vl + 4a 



(27) 



The 1-particle reduced density matrices can then be eas- 
ily obtained by quadrature; for Particle 1, 



/OO 
'iqBp{qA,qB;q'A,QB) 
-00 



(28) 



2Ci - 2C2 



exp[-Ci{q\ + q'2) + 2C2qAq'A], 



where the state is normaHzed to unity and the constants 
Ci and C2 are 



1 + 2a + 3V1 + 4a 



8 + 8Vl + 4a 



and 



a(Vl + 4a - 1) 



(1 + 2a + VI + 4a) 



(29) 



(30) 



From Eq. l(28|) . we can also define the Gaussian charac- 
teristic length a which characterizes the probability dis- 
tribution of a single particle: 



2mijj 



^/T+Aa 



1 + VI -I- 4a 



(31) 



For bipartite Gaussian states, the entanglement of for- 
mation is known exactly [25| . For the ground state of our 
system, the value is 



5 p(^) = - l0g2 1 - W - 

{l-w) 



(32) 



where 



w ■ 



l + 3Vl + 4a + 2[a- (l + 4a)"- (l + 4a)^ 



\ + 2a- VI + 4a 



(33) 
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IV. METHOD 

For the spin system, the entanglement can be cal- 
culated straightforwardly by standard methods; for the 
two-oscillator system, we calculate the von Neumann en- 
tropy S{p'^^^), and hence the entanglement, numerically 
by using two different approaches. 



A. Expansion in a complete set 

We define an orthonormal set of functions, {4>niq)}, 
with support in a region A of configuration space of width 
2a centred at coordinate q: 



q+a 



(34) 



q — a 

A suitable choice is 



1 {q~q)mr . 

- cos( j n IS odd 

a 2a 

1 . .(q — q)mT 

- sm( ) n IS even, (35) 

a 2a 



= if \q- q\ > a 

We then approximate the appropriate post-selected 
density matrix by an expansion in a finite set of the 
functions defined in Eq. (|35l) : as an example, if only 
Alice makes a preliminary measurement to localise her 
particle in the region A, we have: 



N 



Pd 



(36) 



with pmn given by 



rqA+a rqA+a 

Pmn = / dqA / dq'ji 

JqA-a JqA-a 

^:Mp^^\qA:q'A)Mci'A). (37) 



where p^"*-' (g^i; g^) is Eq. (|28|) . We normalize 

p^j^\qA; q'A) by its trace and can then quantify entan- 
glement by calculating the von Neumann entropy from 

this normalized p-^\qA] q'A) ■ Unfortunately the quadra- 
tures in Eq. i|37p must be performed numerically, making 
this approach relatively time-consuming. 

B. Configuration-space grid 

We therefore explored also a direct real-space ap- 
proach, in which we first discretize the configuration 
space into a finite number of measurement "bins", then 
select only those bins that correspond to the regions 
within which Alice's and Bob's respective particles lo- 
calise. For example, consider again the case in which 



only Alice makes a preliminary measurement, if the re- 
gion is qA — o, < qA < qA + a, we divide this space into 
Nb regions with Nb + 1 equally spaced points {qA 's) cov- 
ering the intervals from qA = qA — a to (ja + a- We then 
build the {Nb -I- 1) x {Nb + 1) post-selected 1-particle 

reduced density matrix p'^\qA', q'A) by calculating its el- 
ements Pmn 's from the 1-particle reduced density matrix 
Eq. 

Pmn = P^^HqT.qA) for l<m,n<NB + l. (38) 

As in the other approach, we calculate the von Neumann 
entropy of the normalized p'^\qA] 5^) in order to quan- 
tify the entanglement. 

Note that if on the other hand both parties make a pre- 
liminary measurement, we start from the full 2-particle 
density matrix and apply Bob's restrictions w.r.t his os- 
cillator before we reduce it into the 1-particle density 
matrix for Alice's oscillator. 

We find that results from the two approaches converge 
to the same values as the number of grid points, or the 
number of expansion functions, tend to infinity. Since 
the second (grid-based) approach is much more efficient 
to compute, it has been used for all the results presented 
in this paper. 



V. RESULTS 
A. The spin system 

1. Pure states 

We present results in Fig. [TJ For the spin system 
we consider, entanglement present in the state Eq. ifTSj) 
depends on 9i and 62 with periods of 7r/2, as shown in 
Fig. WLsl). The maximum entanglement is 2 ebits and 
occurs when both pairs of spins are in the Bell state (^1 — 
02 = {2n+ l)7r/4). When 61 = 62 = mr/2, the state 
reduces to all spins either all up or down so completely 
loses any entanglement. 

Now if the restricted region for both Alice and Bob is 
chosen to be the subspace in which the total z-component 
of spin takes the value zero, and we work in the discard- 
ing ensemble so all other states are ehminated, the entan- 
glement properties of the system become very different. 
Fig. [ijb) shows that the entanglement distribution of 
the restricted state has periods of tt instead of 7r/2, and 
the maximum possible entanglement (now 1 ebit since 
the restricted subspaces for both Alice and Bob are two- 
dimensional) is achieved whenever 9i = 62 or 9i + 62 are 
integer multiples of tt so that the restricted state is in 
the Bell state. Note that there is a singularity whenever 
cos 20i cos 262 = 1. 

If we compare Fig. [ija) and Fig. [Hb), it seems that 
in some instances the restricted state has higher entan- 
glement. This is indeed the case as shown in Fig. [UJc), 

where AS' = S{p"j^^) — S{p^^^) is plotted against both 9i 



(a) Unrestricted 



(b) Restricted (discarding ensemble) 



(c) Differences 



FIG. 1: Entanglement (von Neumann entropy S{p^ ')) present in the chosen spin system (a) when the total spins Ms is 
unrestricted, (b) in the discarding ensemble when Ms for each party must be 0. (c) Entanglement difTerences between the two 
cases; AS' = ^(p*^^) - 5(p('^)). 



and 02 . This is an example of the familiar process of en- 
tanglement concentration f2^. [23|. in which some partial 
entanglement is concentrated after chosen local measure- 
ments. Entanglement is not created on average in our 
example because the probability of finding Ms = is not 
100%. Therefore the inequality Eq. (fl7| is not violated. 



2. Mixed states 

Now we perform a similar calculation for the mixed 
state IpTj) , comparing the entanglement (as quantified by 
the negativity) present when the total spins Mg is unre- 
stricted and the entanglement in the discarding ensem- 
ble when Ms for each party must be 0. The results are 
presented in Fig. [21 We choose three values of F for 
comparison; F = 0.3, F = 0.65 and F = 1. 

In Fig. m we plot the variation of the entanglement 
A/'(/5) with F by choosing both 9i and 02 to be 7r/4 (other 
values can be chosen without affecting the qualitative be- 
haviour but entanglement will not vanish at smaller val- 
ues of F). The case when the total Mg is unrestricted po 
is plotted as a solid line, while the case of the correspond- 
ing discarding ensemble po is plotted as a dashed line. 
The entanglement in both cases vanishes at F = 0.25. 
This is similar to what we observed in a previous pa- 
per [2l|: in that case, we showed that the entanglement 
(as quantified by the negativity) of a two-mode Gaussian 
thermal state vanishes at the same temperature regard- 
less of whether the initial state, or the post-selected state 
in the discarding ensemble, is studied. 



B. Two oscillators: the limit of small region sizes 



For the Gaussian system described in Section llll BI the 
entanglement can be evaulated analytically in the limit 
of very small region sizes, following the method described 
in M- 



1. Only Alice's particle restricted 

Suppose only Alice makes a preliminary measurement, 
and determines that her particle is located in a re- 
gion of length 2a centred at coordinate g^, as in §IVI 
<iA — d < qA < qA + OL- In the discarding ensemble, the en- 
tanglement is — h{e) = — [elog2(e) + (l — e) log2(l — e)] 
with 



a(^/l + 4a- 1) 
12(1 2a + VI + 4a)' 



(39) 



Note that this depends only on a and on the parameters 
of the underlying oscillator system; it is independent of 
q. Note also that the entanglement is non-zero for any 
non-zero a, and can be made arbitrarily large (for a given 
small a) by increasing a. 



2. Both particles restricted 

On the other hand, if both parties make measure- 
ments, thereby also restricting Bob's particle to a region 
of length 26 around (Jb , the entanglement is once again 
h{e) but now e becomes 

(H-2a- VTT4^), (40) 

and the concurrence density [HI] is 



rnu! 



1 + 2a - Vl + 4a 



(41) 



Once again, this result depends only on the dimension- 
less coupHng strength a and the fundamental length unit 
(moj)"^/^ of the oscillators; it is again independent of the 
location of the centres of the measurement regions. Later 
we will see that as a and b increase, the entanglement dis- 
tribution gradually changes so that more entanglement 
is located at some parts of configuration space than the 
others. 



7 



(A) F = 0.3 




(a) Unrestricted (b) Restricted (discarding ensemble) (c) Differences 



FIG. 2: Entanglement (negativity Af{p)) present in the mixed state l[2T1l (a) when the total spins Ms is unrestricted, (b) in 
the discarding ensemble when Ma for each party must be 0. (c) Entanglement differences between the two cases; AJ\f = 
M{f)D) — M{po). F determines the "mixedness" of the state; when -F = 1, the state is pure. 



C. Two oscillators: finite region sizes 

1. Only Alice's particle restricted 

For simplicity, we will set to = 1, = 1 and choose 
the Gaussian characteristic length (Eq. (|3T1) ) for an un- 
coupled harmonic system, ct = 1, as our unit of length. 

In this section, we consider the case in which only Al- 
ice makes a preliminary measurement to determine that 
her particle lies within a finite-size region. Suppose that 
the size of this region is 2a and the location of the cen- 
tre of the region is g, the von Neumann entropy S{p'^^^) 



depends on both 2a and q. This is shown in Fig. |4l We 
look at the variation with q first; Fig. [4] along the g-axis 
shows some of the examples. For finite a, the entangle- 
ment is higher if we measure around the centre of the 
wavefunction, where the probability of finding a parti- 
cle is highest, than if we take our measurements further 
away from the the centre of the wavefunction where the 
chance of finding a particle is very low. 

We can understand this variation by examining Al- 
ice's post-selected reduced density matrix in the centre 
of Fig. = 0) and at the edge \q = ±4). At the edge, 
the diagonal elements increase rapidly towards one end; 
the eigenvalues of this density matrix are dominated by 
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FIG. 3: Variation of the entanglement (negativity 7V(p)) with 
-F € [1/16, 1]. is a quantity that determines the mixedness 
of the state as defined by Eq. jHJ. The solid line is for 
the original ensemble po whereas the dashed line is for the 
corresponding discarding ensemble po- Both Q\ and Q2 have 
been set to 7r/4 to produce the plots. 



these terms, resulting in one eigenvalue being close to 
1 and the other eigenvalues being very small. The von 
Neumann entropy will therefore also be small. In con- 
trast, the diagonal elements in the centre case, instead 
of being dominated by a single element at one end, are 
approximately constant. The resulting spread of eigen- 
values leads to a higher von Neumann entropy. 

We would also expect that as the region size ap- 
proaches the total configuration space, the entanglement 
in the discarding ensemble should tend to the entangle- 
ment originally present in the whole system; this is shown 
in the upper part of Fig. [U where the entanglement rises 
with a until it saturates to the peak value of magnitude 
^ 0.702 given by Eq. Roughly speaking, 

this saturation occurs once the region has expanded to 
include a significant portion of the central part of the 
harmonic oscillator wavefunction. 

We have already seen that in the Hmit of small a the 
entanglement becomes independent of position. In fact, 
even for finite a the entanglement is distributed very dif- 
ferently from the probability distribution of Alice's par- 
ticle. This is shown in the lower part of Fig. IH where 
the coloured curves show the entanglement (scaled to a 
common maximum value) as a function of q for different 
widths 2a; for comparison, the black dashed plot shows 
the Gaussian one-particle probability distribution with 
standard deviation a given by Eq. (|3T|) . Note that the 
width of the entanglement plot varies non-monotonically 
with a: the entanglement is constant in the limits of small 
and large a, and has a minimum width around 2a = 2 
(for a = 6). Note also that S{fl^^'^) is very small but is 
non-zero even for small a, as expected from Eq. (|39|) . 

For comparison, we also present in Fig. [DJb) results 
for a much weaker coupling, a — 0.06 compared with 
(2 = 6: for weak coupling, the entanglement has smaller 
peak values (= 0.00859 in this case) and its spread is 
narrower, but the qualitative features are similar in both 
cases. 



2. Both particles restricted: entanglement distributions 

Next we consider the case where both Alice and Bob 
make preliminary measurements, but not necessarily in 
the same way. 

We start by considering two different cases; the first 
(Case 1) is that both parties' preliminary measurements 
restrict their particles to regions with identical widths 
and centres {a = b and qa = Qb), whereas in the sec- 
ond case (Case 2) the region widths are the same but 
the centre of Bob's region is always fixed around the cen- 
tre of the wavefunction (a = 6, qb = 0). The results, 
for a = 6, are shown together with the previous result 
(Case 3; only Alice makes a preliminary measurement, 
as shown in Fig. IDJa)) for comparison in Fig. [5j The 
entanglement in the discarding ensemble of Case 3 is the 
highest out of the three cases; this is as expected, since 
the entanglement can only reduce under the additional 
(local) measurements made by Bob. When the width 2a 
is small, the entanglement of Case 1 is higher than of 
Case 2. However, as 2a increases. Case 2 converges more 
rapidly to Case 3 so that its entanglement is now higher 
than that of Case 1, until 2a becomes so large that the 
differences between all three cases disappear. 



3. Both particles restricted: classical correlations 

We now compare the entanglement distributions to the 
classical correlations between the particles. Suppose that 
Alice and Bob localise their respective particles to re- 
gions with the same widths but different centres; the en- 
tanglement in the discarding ensemble will depend on 
both qA and qs- We shall compare the entanglement 
distribution with the 2-particle probability distribution 
P{qA & Anqs € B), and the conditional probability dis- 
tribution for Bob's particle given the position of Alice's 
particle, P{qB & B \ qA & A). 

The two-particle probability is 



MA+a rqB+a 
P{qA eAnqeeB) = / dqA / 



dqi 



p{qA,qB;qA,qB), (42) 

and in the limit of small a, 6 we have 

P{qA e AnqB e B) = Aabp{qA, qB;qA, m)- (43) 
The conditional probability is 

P{qB eB\qAeA)^ P(^Z^) ' ^^^^ 

where P{qA G A) is the 1-particle probability. In the 
limit of small a, b this becomes 

ut ^ K \ ^ A\ nhP{^^i^B]qA,qB) ..p.^ 

P{qB (^B\qA(^A)^ 2b rp-— — ^— . 45 

P^^'{lA,qA) 
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(a) a = 6 (b) « = 0.06 

FIG. 4: (Color online) Top: Variation of the entanglement 5'(p''^') with both the width 2a and the centre q of the preliminary- 
measurement region. Bottom: S{p^^^) plotted against q for different widths, re-scaled such that S{p^^'^) has the same peak 
value at g = 0. A plot (the black dashed line) of the corresponding Gaussian probability distribution for Alice's particle, with 
a standard deviation a determined by the coupling strength a, is shown for comparison. The different plots correspond to two 
different coupling strengths, (a) a = 6 v.s. (b) a = 0.06. The number of bins Nb used in the calculation was 200 in both cases. 




(a) 2a = 0.5 (b) Larger 2a's 



FIG. 5: (Color online) Comparison of the two different cases of preliminary measurements done by both parties together with 
the case that only one party makes a preliminary measurement. The entanglement S{p^^^) is plotted against the centre q of 
the preliminary-measurement region with width 2a. (a) For 2a = 0.5. (b) For other larger values of 2a. Red long-dashed line 
(Case 1): Both parties' preliminary measurements localise their particles in regions with identical widths and centres {a = b 
and cjA = qs). Blue thick short-dashed line (Case 2): The widths of the regions are the same but one centre is always fixed 
around the centre of the wavefunction while there is no restriction on the other centre (a = b, qs — 0). Black thin solid line 
(Case 3): Only one party makes a preliminary measurement. In all three cases, the number of bins used in the calculation is 
Nb = 100 and a = 6. 
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In each case the small-a, b limit can be easily evaluated: 
we find 

P[qA,qB;qA, qs) = C,2 exp( — ^ ) 



2<jI 



with 'classical' standard deviations 



2ai 



CTI = \/2, 



c 



and 

p{qA,qB;qA,qB) 



p^'^KqA.qA) 



C3exp(-^ 



qAqs 



(46) 



(47) 



ft) 



with 



0-2 

fl2 



^ (1 + Vl + 4a)(l + 2Q + Vl + 4a) ^i ^^^^ 



4a^ 



1 + Vl + 4a 
1 



Vl + 4a- 1' 



where (^2 and C3 are normalization constants. 

For finite a and b we capture the shape of the distri- 
butions by fitting the numerically calculated values of 
P{qA € Ar^qB & B), and P{qB £ B \ qA & A) using 
the same expressions Eq. (|46| and Eq. l(48|) . thereby ex- 
tracting numerical values for cr^, cri,2 and a 12- We also 
use the function Eq. (|46l) to fit the entanglement dis- 
tribution, thereby obtaining two further parameters cr^ 
which quantify the extent of the entanglement distribu- 
tion along its principal axes. 





(a) 2a = 0.5 



(b) 2a = 4 



FIG. 6: Dependence of the entanglement 5'(p^'*') on the loca- 
tions of the centres of the preliminary-measurement regions 
QA and qB- (a) the width 2a of the regions is 0.5. (b) 2a = 4. 
In both cases, Nb = 100 and a = 6. 

As before, we take a = 6. In Fig. [6l we show two 
cases of entanglement distributions for different widths 
(2a ~ 0.5 and 2a — 4) of the preliminary-measurement 
regions. We see that the entanglement distribution with 
larger 2a is more symmetric. The corresponding joint 
probability distributions and conditional probability dis- 
tributions are shown respectively in Fig. [7| and Fig. [8l 





(a) 2a = 0.5 



(b) 2a = 4 



FIG. 7: The dependence of the classical joint probability 
P{qA e An qb € B) on cjA and qb. (a) 2a = 0.5. (b) 
2a = 4. In both cases, a = 6. 




(a) 2a = 0.5 



(b) 2a = 4 



FIG. 8: Dependence of the conditional probability P{qB £ 
B \ QA € A) on QA and qb- (a) 2a — 0.5. (b) 2a = 4. In both 
cases, a = 6. 



(Note that the figures show different range of (ja and 
(jB-) The classical probability distributions P{qA H qs) 
are more localized and symmetric in space than the en- 
tanglement distributions. 



a = 6 












0-2 


0"12 


2a ^ 


00 


00 


1.41 


0.632 


0.866 


0.577 


0.500 


2a = 0.5 


10.4 


2.29 


1.43 


0.665 


0.937 


0.603 


0.531 


2a = 4 


3.44 


2.10 


2.37 


2.00 


11.0 


1.53 


2.64 



TABLE I: Table of a values for a = 6. 



In the limit of very small a, ^(pf^)) is constant ev- 
erywhere (Eq. (|40l) ) so (J^ and cr? must diverge; the 
results in Table U show that ct^ diverges more quickly as 
a reduces, while the two parameters become compara- 
ble for large a as the entanglement distribution becomes 
more symmetric. Indeed, the distributions of the entan- 
glement and the classical correlations become more alike 
as 2a increases, because both distributions are flat out 
to a distance a either side of the wavefunction's central 
peak. 

We can also study the effect of varying the coupHng 
strength a for a fixed (small) 2a. We plot and cr? 
against a with 2a — 0.5 in Fig. [9K whereas a^, a^, (Ti, 
(72 and (712 in Fig. [Sb- The entanglement distribution is 
the most asymmetrical and as a increases, the difference 
between and a'^ widens. Of the quantities deter- 



mining the classical probability distribution, cr^ remains 
constant with increasing a, but cr^ gradually decreases. 
These trends arise because the two particles tend to move 
together when the spring joining them becomes strong. 
Therefore, as a increases, the white rod in Fig. [8] rotates 
about the centre of the square from the line qb = to- 
wards the diagonal qa — Qb- cti is always the largest out 
of the three parameters for the conditional probability 
distribution. For weak a, ai2 is larger than (72 but as 
a becomes larger, at some point the two plots intercept 
and tTi2 is no longer larger than a2- 

How in the limit of very small a these quantities (Eq. 
(|47l) and Eq. (1491) ) vary with a is shown in Fig. [91:. We 
see that the behaviour of these quantities do not change 
much, compared with the previous results when 2a = 0.5, 
apart from that the interception points happen at smaller 
a. Note that a^. diverge as a ^ 0, so these parameters 
are not shown. 



Willi h 




5 ID 15 20 Z5 30" 



(a) crj; 2a = 0.5 



Width 




5 10 15 2D 25 3D " 



(b) 2a = 0.5 




5 10 15 20 26 IS " 



(c) Small-a limit 

FIG. 9: (Color online) Plots of (T+, (t_, cti, ai2 and (T2 against 
a. In the plot legend, Q stands for the 'quantum' entangle- 
ment distribution and C for the 'classical' probability distri- 
bution, (a) and (b): Numerical results: 2a is chosen to be 0.5 
for all the cases, (c) Analytical results: in the limit of very 
small a. 
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VI. CONCLUSIONS 

We have presented a novel thought experiment that 
gives an approach to determining the location in configu- 
ration space of the entanglement between two systems. It 
involves choosing a region of the two-party configuration 
space and making a projective measurement with only 
enough resolution to determine whether or not the sys- 
tem resides in this region, then characterizing the entan- 
glement remaining in the corresponding sub-ensemble. 
Our approach is particularly simple to implement for 
pure states, since in this case the sub-ensemble in which 
the system is definitely located in the required region af- 
ter the measurement is also a pure state, and hence its 
entanglement can be simply characterized by the entropy 
of the reduced density operators. 

We have given examples of the application of our 
method to states of a simple spin system, where Alice 
and Bob share two pairs of spin- 1/2 particles, and also of 
a continuous-variable system in which they share a pair 
of coupled harmonic oscillators. 

The first case shows how the amount of entanglement 
located in the chosen region (in this case the ~ 
manifold) varies as the characteristics of the states shared 
by Alice and Bob are altered. We presented results for 
both pure and mixed states, and show how entanglement 
is affected by the "mixedness" F both qualitativly and 
quantitativly. Specifically, we show that the states which 
are entangled from the global point of view are also en- 
tangled by our local measures [2l|, i.e. global entangle- 
ment of the initial state vanishes at the same point as the 
entanglement remaining in the discarding ensemble after 
the preliminary measurement to locate the system in a 
chosen subspace. 

For the second case we have presented results as a func- 
tion of the strength of the coupling between the oscilla- 
tors, as well as of the size and location of the prelim- 
inary measurement regions. In all cases the remaining 
entanglement saturates to the total entanglement of the 
system as the measured regions become large. For small 
measured regions the entanglement tends to zero, but for 
a fixed region size the configuration-space location can be 
varied in order to give a variable-resolution map of the 
entanglement distribution. We find that the distribution 
of the entanglement is qualitatively different from the 
classical correlations between the particles, being con- 
siderably more extended in configuration space than the 
joint probability density and becoming more and more 
diffuse as the size of the regions decreases. 

Our approach suffers from the disadvantage that there 
is no sum rule on the entanglements in the discarding 
ensemble: the sum of the entanglements from all the 
sub-regions defined by a given decomposition of configu- 
ration space does not yield the full entanglement of the 
system. Instead, the entanglements from the sub-regions 
satisfy the inequality in Eq. iflTl) . It would be interesting 
to understand in more detail the relationship between 
the restricted entanglements (as defined in this paper) 
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and the full entanglement of the system, and also to ex- 
tend the calculations reported here to projective mea- 
surements made in other bases, to POVMs, and to mixed 
states. 
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